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Abstract 

We study interacting particle systems on the real line which generalize 
the Hammersley process [D. Aldous and P. Diaconis, Prob. Theory Relat. 
Fields 103, 199-213 (1995)]. Particles jump to the right to a randomly cho- 
sen point between their previous position and that of the forward neighbor, 
at a rate which may depend on the distance to the neighbor. A class of 
models is identified for which the invariant particle distribution is Poisson. 
The bulk of the paper is devoted to a model where the jump rate is constant 
and the jump length is a random fraction r of the distance to the forward 
neighbor, drawn from a probability density </>(r) on the unit interval. This 
is a special case of the random average process of Ferrari and Pontes [El. J. 
Prob. 3, Paper no. 6 (1998)]. The discrete time version of the model has 
been considered previously in the context of force propagation in granular 
media [S.N. Coppersmith et al, Phys. Rev. E 53, 4673 (1996)]. We show 
that the stationary two-point function of particle spacings factorizes for any 
choice of (j){r). Under the assumption that this implies pairwise indepen- 
dence, the invariant density of interparticle spacings for the case of uniform 



(j){r) is found to be a gamma distribution with parameter z^, where = 1/2, 
1 and 2 for continuous time, backward sequential and discrete time dynam- 
ics, respectively. A heuristic derivation of a nonlinear diffusion equation is 
presented, and the tracer diffusion coefficient is computed for arbitrary 0(r) 
and different types of dynamics. 

KEY WORDS: Interacting particle systems; random average process; invariant 
product measures; discrete time dynamics; hydrodynamic limit; single file diffu- 
sion; granular packings. 

1 Introduction and outline 

In this paper we are concerned with systems of interacting particles moving on the 
real line. The models of interest can be described as follows: Let Xj G M denote 
the position of the i-th particle. In an elementary move particle i jumps to the 
right to a position Xj + 6i between Xi and Xj+i > Xj. In the absence of a lattice 
spacing, there are two natural ways of setting the scale for the jump distance 5,: It 
can be imposed externally through the choice of a fixed probability density fi{Si), 
in which case moves with 6i > Xj+i — Xi have to be rejected, or the scale can be 
set by the gap or "headway" 

Ui = Xi+i-Xi (1.1) 

in front of particle i by letting fi depend on the configuration U = {ujjjgz as 

m\U)=ur'<f){6,/ui), (1.2) 

where 0(r) is a probability density with support on the unit interval. Equation 
( p..2|) implies that the jump length 6i is a random fraction r of the headway Ui. 
The rate for the move is a function 'y{ui) of the headway. The moves are executed 
in continuous time (in which case each particle is equipped with an exponential 
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clock) or in discrete time; in the latter case the particle positions are updated either 
in parallel, or sequentially by going through the system against the direction of 
particle motion. A model is defined by specifying the functions fi{S) and 7('u) as 
well as the type of dynamics (continuous time, parallel or sequential). 

Two equivalent representations of the dynamics will prove to be useful. In 
terms of the headway variables Ui the particle configuration may be visualized as 
a system of sticks located at the sites i of the integer lattice, Ui being the length 
of stick i. In an elementary move a fraction 6i of stick i is broken off and added 
to stick i — 1 0, Alternatively, the particle positions Xj(t) can be taken to 
define the height of a one- dimensional interface above the point i. The asymmetric 
particle motion translates into a growth process, and the fact that particles cannot 
pass each other implies that the interface is a monotonically increasing staircase 
{xi^i — Xi > 0) at all times. We will refer to these two viewpoints as the stick 
representation and the interface representation, respectively. 

For continuous time dynamics, a jump length distribution of the type (|1.2| ) with 
uniform, and '~f{u) = u the model reduces to the Hammersley process discussed 
in [|T|. In this case the invariant distribution of particle positions is Poisson. Here 
we are interested in obtaining similar results for other choices of fi and 7, and 
other types of dynamics. Our motivation is mainly conceptual: While a wealth of 
results [E^, 521 ITR EBI, B2| are available for particle systems on the integer lattice 



such as the asymmetric simple exclusion process little is known analytically 
for the case of continuous particle positions, although motion on the real line 
appears naturally e.g. in applications to highway traffic |]TB|, |TB|. 

An important simplifying feature of the asymmetric exclusion process is the 
existence of stationary product measures. Here the analogous desirable property 
is the product form 

V{U)=\[P{u,) (1.3) 
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for the stationary probability of a configuration U of particle headways. Therefore 
a primary goal will be to find nontrivial examples of asymmetric particle systems 
on R for which ( |1.3|) holds. 

We provide an outline of the paper. In the next section we explore the condi- 
tions for a Poisson distribution of particle positions (corresponding to an exponen- 
tial distribution of interparticle spacings in ( |1.3|) ) to be invariant for continuous 
time dynamics. Our strategy is to consider a finite number of particles mov- 
ing on a ring of length L, and to demand that the stationary measure gives the 
same weight to all allowed configuration; this then implies a Poisson measure for 
A^, L — > oo at fixed density p = N/L. Provided the jump rate 7 is independent of 
the headway, we find that the Poisson measure is invariant for arbitrary externally 
imposed (i.e., configuration and particle independent) jump length distributions 
f{u). On the other hand, if the jump length is scaled to the headway as in ( p..2| ), 
the Poisson measure is stationary only for a one-parameter family of power law 
functions and 7, which have been identified previously in the context of (sym- 



metric) stick models JTO 



Sections 3 and 4, which constitute the main part of the paper, are devoted to 
models with constant jump rate, 7 = 1 independent of the headway, and jump 



length distributions of the type ( |1.2| ). In the interface representation these belong 
to the class of random average processes (RAP) studied by Ferrari and Pontes |12[| : 
The particle position x'^ after the move is an average 

x[ = rxi + {1 - r)xi+i (1.4) 

of the previous positions Xi, Xj+i, with a random weight r G [0, 1] drawn from 
the probability density 0(r). We therefore refer to these models as Asymmetric 
Random Average Processes (ARAP). Discrete time ARAP's have been introduced 
previously to model force fluctuations in random bead packs ||2^, 0, ||. In that 
context the headway Ui{t) represents the (scaled) force supported by bead i at 



depth t below the surface of a two-dimensional packing (see Section 3.2.1). 

In Section 3.1.1 we show, for the case of continuous time dynamics, that the 
two-point correlation function of particle headways (uiUj) factorizes in the station- 
ary state for any choice of 0(r), and obtain the expression 



of 0(r) and the particle density p. Similar results for the discrete time models are 
derived in Section 3.2. 

More detailed information about the stationary headway distribution can be 
obtained when 0(r) is the uniform distribution on [0, 1]. Assuming that the fac- 
torization property of the two-point function implies pairwise independence of the 
Ui, we derive and solve stationarity conditions for their moments, which show that 
the invariant density of headways (normalized to (ui) = 1) takes the form of a 
gamma distribution, 



where the parameter u depends on the dynamics: For continuous time dynamics 
u = 1/2, while sequential and parallel dynamics yield u = 1 and 2, respectively. 
The result for parallel dynamics has been previously derived by Coppersmith et al. 
0, who also gave an explicit proof of the factorization property (p..3|). Equation 
( |1.7|) implies hunching of particles (enhanced density fluctuations compared to the 
Poisson measure) for continuous time dynamics (z/ = 1/2) and antibunching for 
parallel dynamics (z/ = 2). The associated nontrivial particle-particle correlations 
are explicitly computed in Section 3.3. 

Based on numerical simulations, we conjecture that the stationary single par- 
ticle headway distribution is exactly given by ( |1.7| ) for all three types of dynam- 




(1.5) 



for the stationary variance of headways in terms of the moments 






(1.7) 
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ics. For continuous time dynamics and a finite number of particles on a ring tlie 
assumption of an invariant product measure is examined in Section 3.1.2. Surpris- 
ingly, we find that the product measure is not invariant for the ARAP, although 
it is invariant for a related symmetric stick model. This conclusion agrees with 
recent results for the infinite system obtained by Rajesh and Majumdar [^. 

Section 4 is devoted to the large scale, long time behavior of the ARAP. We 
derive a hydrodynamic equation of singular diffusion type, and compute the tracer 
diffusion coefficient using a Langevin approach. Since these results depend only 
on the stationary two-point function of headways, they are valid for any choice of 
the jump length distribution 0(r). Finally, some conclusions and open questions 
are formulated in Section 5. 



2 Models with invariant Poisson measures 

2.1 Constant invariant measure on the ring 

In this section we want to identify continuous time dynamics which leave a Poisson 
distribution of particle positions invariant. For this purpose we first consider 
particles moving in continuous time on a ring of length L, with density p = N/L. 
Allowed headway configurations then satisfy the constraint 

N 

^Ui = L (2.1) 

i=l 



and the product measure (O) is required to hold on the set of configurations 



defined by (|2.1j ). For an exponential distribution P{u) ~ e this implies that all 

allowed headway configurations carry the same weig ht Vl{N,L)-^, where 

tN-i 

denotes the volume of the set, i.e. the invariant measure is constant on allowed 
configurations. It is straightforward to check that this implies Poisson measure in 
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the hmit A^, L ^ oo at fixed density p. For example, the distribution of a single 
headway on the ring is given by 

i-^.LW = "^^^"^^;;^p^ --pe-^", iV,L^oo (2.3) 

while the joint distribution of the headways of two neighboring particles is 

Pn,l{u^, = " ~ "'^^^ - p2^-PK+n.,o AT, L ^ oo. (2.4) 

A similar argument can be carried out for the probability distribution of the par- 
ticle positions on the ring. 

Invariance of the constant measure requires the total transition rates for going 
into and out of any configuration to balance. This yields the condition 

V/ ~ dwMw\U^'\w))^{u, + w) = y2 dwUw\l{)^iu,) (2.5) 
i=i i=i 

for any configuration U, with the configuration U^'^\w) = {u^j\w)}j(zz defined 

through 

f 

Ui + w : j = i 

-w : j = i -1 (2.6) 



u)'{w) = < 



Uj : else 



and periodic boundary conditions implied in the summation over i. Note the upper 
integration limits, which ensure that particles cannot pass each other {6i < Ui). 
Two examples of dynamics which satisfy (|2.5| ) will be given in the following. 



2.2 Configuration-independent jump length distributions 



If the jump rate 7 is independent of headway, the invariance condition (|2.5|) is 
seen to hold for any jump length distribution f{w) which is independent of the 
configuration and of the particle label i. The stationary speed v of particles at 
density p is then computed from 

v = -ip due-P^ dwwf{w), (2.7) 
Jo Jo 



and the current follows from j{p) = pv{p). For example, for jump lengths chosen 
uniformly in the unit interval one finds 

j(p) = ^[l-(l + p)e-1. (2.8) 

It should be noted that in general the Poisson distribution is not the unique 
invariant measure. For example, if fiw) = for w less than some minimum 
jump length a, then all configurations with Ui < a for all i are trivially invariant. 
Numerical simulations indicate, however, that such "absorbing" states are typically 
not reached, even if the system is started very close to them. If f{w) = 6{w — 
1) and the particles are started on the integer lattice, the model reduces to the 
asymmetric exclusion process, which has a geometric (rather than exponential) 
headway distribution. 



2.3 Scale-invariant models 



When the scale of the jumps is set by the headways, inserting ( |1.2| ) into (2!^) and 



requiring the terms on both sides to cancel pairwise yields the following integral 
equation connecting the functions (f) and 7, 

■ , , .(l)(w/(u' + w)) , , ,^ 
u' + w 

which should be true for all u, u'. Taking the derivative with respect to u this 
becomes a differential equation for 7, 

c[7 ^ -f{v)(f){u/v) 
du V 

with V = u + u' > u. Setting in particular v = u we see that 7 has to be a power 
law function, 

7(w) = 7o«"-\ (2.11) 



8 



where 70 > is a constant and a = 1 + 0(1). Using ( |2.10|) the jump length 
distribution is then found to be also a power law, 

cj){v) = (a - (2.12) 

Normalizability of requires a > 1. 



Equations (|2. 111 , 12^1^) define a one-parameter family of models for which the 
Poisson distribution of positions is invariant for an arbitrary number of particles 
A^, the Hammersley process being given by a = 2. The corresponding symmetric 
stick models, in which the broken-off piece is distributed with equal probability to 
the left or right neighbor, were considered by Feng et al. |10|- Since 7 is a power 



law, these models are scale invariant in the sense that the average particle spacing 
(wj) = 1/p is the only length scale in the problem. Therefore also the stationary 
particle current j is a power law function of the density. To compute it, we note 
that the average particle speed is given by 

/■oo pi 

V = {j{ui)5i) = p due-'^j{u)u dvv^{v) = -fo{l-l/a)T{a + l)p-'^ (2.13) 
Jo Jo 

and therefore 

j{p) =pv = 7o(l - l/a)T{a + l)pi-". (2.14) 



3 Asymmetric Random Average Processes 

The asymmetric random average process is a scale-invariant model characterized 
by a jump length distribution of type (p..2[), and a constant jump rate 7 = 70 = 1. 
The discussion is phrased most naturally in the stick representation, and begins 
with the continuous time models. 



3.1 Continuous time dynamics 
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3.1.1 Stationary headway correlations 

Consider first the time evolution of the second moment {uf). In a small time 
interval At two processes affecting Ui may occur: A random fraction 5, of Ui may 
be lost to z — 1, and a random fraction of Ui+i may be gained from i + Both 
processes occur with probability At. Thus 

{ul){t + At) = At[((«, - 5,)') + {{u, + k+if)] + (1 - 2At)(«^)(t). (3.1) 

Stationarity then implies 

- 2{5,u,) + {6f) + 2(5,+iM,) + (52^1) = 0. (3.2) 

Since 5j = TjUj where rj is an independent random variable with mean fii and 
second moment fi2, we have that {6iUi) = fii{u^), {5f) = = fJ'2{uf) and 

{6i+iUi) = fii{uiUi+i). Thus (|3^ ) becomes 

{(Xl - /i2)(Mi) = Hl{UiUi+i). (3.3) 

Similarly for the general two-point function = {uiUi^k) we obtain the station- 
arity condition 

/xi(Cfc+i + Ck-i - 2Ck) = /i2Co('5fc,i + h-i - 2(5fc,o), (3.4) 

where translational invariance and symmetry (C^ = C_fc) of the correlations has 
been used. Solving eg. ( |3.4| ) starting from k = one finds 

Cfc = [l-(/X2^i)(l-M]Co. (3.5) 

Imposing the boundary condition limfc_,oo Ck = {ui)"^ = l/p^ for an infinite system 
of density p, eq.( p.5|) then shows that the two-point function factorizes for any 



A; > 1 and the variance of headways is given by ( |1 . 5|) . 
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3.1.2 Stationary headway distribution for uniform 0(r) 



We now specialize to the case when the distribution of scaled jump lengths 0(r) 
is uniform in [0, 1], and assume that the factorization property which was verified 
above for the two-point function implies the pairwise independence of the Ui. Then 
the stationarity condition for the n-th moment 



yields (the index i of Ui is now dropped) 

n _ 
k=0 

which can be rewritten as a recursion relation, 



(3.6) 



(3.7) 



(3.^ 



k=l 



Evaluating this expression for n = 1, 5 we find that the relation 



n(2A:-l) 



,fe=i 



{uY 



(3.9) 



appears to hold, which is characteristic of the gamma distribution ( |1.7|) with pa- 
rameter u = 1/2. 

To prove it, we first insert ( ^.9] ) into ( |3.(j| ), and obtain 



2n\ _n + l 
n I n — 



1^ 1 f2k\ 

1 ^ A; + 1 V A; / 

k=i ^ ^ 



1 f2k\ f2{n~k) 
n ~ k 



(3.10) 



This can be verified using the binomial expansion 



k=l ^ ' 



(3.11) 



Integrating with respect to x we also have 



1 



;i - 4x) 



1/2 



1 X 

-4 + 2 



+|:^(t-i>-- (3.12) 
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Since the product of the left hand sides is a constant, all coefficients of x"^ with 
m > in the series obtained by multiplying ( |3.11| ) and ( |3.12| ) must vanish. After 



rearranging terms this is seen to imply ( p.lO|) . 



In fact the relation ( |3.9| ) was first guessed on the basis of numerical simulations. 
Rather accurate numerical estimates for the stationary moments of Ui can be 
obtained by starting from an ordered initial condition [ui = 1 for all i) and fitting 
the finite time data to the form 

(M")(t) = A„ + 5„ri/2 (3.13) 

which is suggested by the fiuctuation theory of Section 4.2 (see eg. ( [4.21| )). The re- 
sults shown in Table I strongly indicate that the stationary single particle headway 
distribution is exactly given by the z/ = 1/2 gamma distribution. 

To test the assumption of an invariant product measure underlying the deriva- 



tion of (|3.8| ), we proceed as in Section 2 and consider a finite number of particles 



on a ring. The condition for the product measure ( |1.3| ), restricted to the set (|2.1| 
of allowed configurations, to be invariant now reads 



f /"'"' + = f = A,. (3,14) 

Ui + w P(iii)F(t,i_,) 

Inserting the gamma distribution with parameter u = 1/2 (eg. ( |1.7|) ) and noting 
that 

dw{u + w)-^/\v - wy'/' = (3.15) 

U + V 

the condition (|3.14|) becomes 







with periodic boundary conditions, U-i = un- Eguation (|3.16|) is satisfied for 



N = 2, but not for general A^. We conclude that the product measure (|1.3|) is 
not invariant for A^ different from 2. It is however the exact invariant measure for 
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the symmetric stick process obtained by transferring the piece broken off stick i 
to i — 1 or i + 1 with equal probabihty. Indeed, in that case the left hand side of 
( p.l6| ) becomes 



While these arguments are restricted to finite systems, the conclusions agree with 
calculations carried out for the infinite system by Rajesh and Majumdar ||26|| . 
Specifically, they show that the product measure ansatz for the continuous time 
ARAP breaks down at the level of three-point correlations, but is exact for the 
symmetric stick model. 

3.2 Discrete time dynamics 
3.2.1 Parallel update 

A discrete time version of the ARAP is obtained by writing 



where 6j = rjUj with independent random numbers vj distributed according to 
the density 0(r). This is closely related to a model introduced by Coppersmith, 
Liu, Majumdar, Narayan and Witten for the description of force fluctuations in 
bead packs [0. To see the connection, let W{i,t) denote the weight supported 
by bead i in the t-th layer below the (free) surface of the packing. The key 
assumption of the model is that the beads are arranged on a regular lattice, and 
that each bead transfers its weight to exactly M beads in the layer below. The 
fraction qij{t) G [0, 1] of the weight of bead i in layer t which is transferred to 
bead j in layer t + 1 defines a matrix with random entries subject to the constraint 
J2j lijit) = 1- Assigning unit mass to each bead, the weights evolve according to 



N 




(3.17) 



Ui{t + I) = u^{t) - Si{t) + Si+i{t) 



(3.18) 




(3.19) 
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For large t all weights increase linearly with t, which suggests to introduce normal- 
ized variables U{i,t) = W{i,t)/t. Specializing to a two-dimensional lattice where 
the beads are labeled such that bead i is connected to beads i and i + 1 in the layer 
below, we see that for t oo the evolution of the U{i, t) reduces to ( |3.18| ) with 
the identification qu = 1 — rj and gj+ij = rj+i. In the context of beak packs qu and 
qi+ii should have the same distribution, and hence strict equivalence between the 
two models holds only when 0(r) is symmetric around r = 1/2. 

Let us first show that the stationary two-point headway correlations factorize 
for any 0(r). Proceeding as above in Section 3.1.1, we obtain the stationarity 
condition 

{fii - + Ck-i - 2Ck) = (/i2 - /i?)Co(4,i + 4,-1 - 24,o), (3.20) 

with the solution 

Ck = [l- (/i2 - f^D/ifii - /i?)(l - h,o)]Co. (3.21) 

As in the continuous time case this implies factorization for A; > 1 in the infinite 
system, with the stationary variance of headways given by 

in') - {nf = ,f ^ . (3.22) 

For the case of a uniform distribution 0(r), Coppersmith et al. (see also 
25, |2^) have shown explicitly that the stationary measure takes the product form 
(p..3|), with the headway distribution P{u) given by the gamma distribution ( |1.7| ) 
with V = 2. The latter is easily derived along the lines of Section 3.1.2. Under 
the assumption of pairwise independence, the stationarity condition for general 
moments (m^) now reads 

(^") = 7 — ^r— frt?) (^"-')(«')- (3.23) 

(n- l)(n + 2) ^ + 1/ 
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A straightforward computation shows that this is solved by the expression 

(m") = 2-"(n + (3.24) 
for the moments of the gamma distribution ( p..?] ) with parameter u = 2. 

3.2.2 Ordered sequential update 

In the context of traffic modeling ^ it has been found useful to implement 
a different kind of discrete time dynamics, in which the particles are moved one 
by one, in the order of their positions in the system. This ordered sequential 
update can proceed either in the direction of particle motion (forward update) or 
against it (backward update). For the ARAP it is easy to see that the forward 
update is equivalent to the parallel dynamics discussed in Section 3.2.1, however 
the backward update is not. 

In the stick representation, backward sequential update implies that stick i 
first receives a random fraction of stick i + 1, placing it in an intermediate state 
of length -u^, and subsequently transfers a random fraction 6'^ of to stick i — 1. 
It is important to note that, at the time of transfer of mass to stick i, stick i + 1 
has already received mass from i + 2 and thus the amount transferred from i + 1 
to i is a random fraction of > Uj+i. The dynamics therefore proceeds in two 
steps, 

u[it) = u,it) + 5:^,it) (3.25) 

M,(t+i) = M:(t)-5,:(t), (3.26) 

where Sj is a random fraction of u'j. Taking the average of both sides of ( p.25|) or 
( 3.26 ) the stationary mean of m- is seen to be 



in') = = -jT^y (3.27) 

1-/^1 p(l-/ii) 

Equation (|3.26|) implies the relation 

Ck = [(1 - fXiY + (/i2 - f^l)h,o]C', (3.28) 
15 



between the stationary two-point functions Ck of Ui and of u[. Using (|3.25|) it 



is easy to show that the stationarity condition for is identical to the condition 
( p. 201 ) obtained in the case of parallel update. Therefore also factorizes in 
the infinite system, and through (|3.28| ) this property carries over to Ck- For the 
stationary variance of the backward sequential update model we find the expression 

{u') - {u? = (3.29) 

p^(l -/ii)(/ii -/ia) 

Turning to the stationary headway probability distribution P{u), we again 
assume pairwise independence and note the functional equation 

P(u) = [ dr r~V(l - r)P'{u/r) (3.30) 
Jo 

relating P{u) to the distribution P'{u') of the intermediate state headway. For 
uniform 0(r) the stationarity condition for the n-th moment of u'^ then reads 

(KD = {{u..+5[^,r) = Q (3.31) 

k=0 



Using the relation ((«')") = {n + l){u^) obtained from (|3.30| ) this reduces to 



(n") = ^5^(^)(n^)(n"-'=), (3.32) 

k=0 

which is solved by setting (m") = n\{u)"'. We conclude that P{u) is an exponential 
distribution (a gamma distribution ( |1.7|) with i> = 1). This is confirmed by the 
numerical data shown in Table I. 

From ( |3.30 ) the distribution of the intermediate state headway is found to be 



a z/ = 2 gamma distribution with mean 2{u) = 2/p, 

P\u) = p^ue-P^. (3.33) 

Given the equivalence between the intermediate state headway and the headway 
for parallel update which we found on the level of the two-point function, it is no 
surprise that ( |3.33|) is identical, up to a scale factor, to the headway distribution 
P{u) for parallel dynamics. 
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3.3 Particle-particle correlations 

In this section we illustrate how the product measure (|1.3| ) with the headway 
distribution (|1 . 7| ) translates into nontrivial particle-particle correlations when v ^ 
1. For example, the probability density g{x) for finding a particle at x, conditioned 
on having a particle at the origin, can be written as 

oo 

g{x) = Y,Pn{x), (3.34) 

n=l 

where Pn{x) is the probability density for the n-th particle to be at x when the 
0-th is at the origin or, equivalently, the probability that Y^^=o = ^- The P„ 
are obtained iteratively from Pi{x) = P{x) through the convolution 

Pn{x) = / dy Pn-i{y)P{x - y). (3.35) 

Inserting the gamma distributions ( |1.7| ) with parameters z/ = 1/2 and = 2, one 
finds that 

Pn{x) = p(r(n/2)2"/2)-i(pa;)"/'-^e-''"/2 (3.35) 
for the continuous time case, and 

Pn{x) = ^^;^(P^)'"~'e"'''^ (3-37) 

for parallel dynamics. 

In the parallel case the evaluation of the sum ( |3.34| ) is straightforward, and 
yields the expression 

g{x) = p(l - e-'n (3.38) 

for the correlation function, which explicitly displays the tendency of particles to 
avoid each other at distances short compared to 1/p. 

To compute ( |3.34| ) with the Pn given by ( p.36| ), it is useful to write g as the 
sum of two contributions (/even and godd from even and odd n, respectively. One 
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finds that 5'even(a;) = p/2 independent of x, wliile tlie odd part can be brouglit into 
tlie form 

^ m=l ^ 

To sum tlie series we write (m — 1)! = dz z"^~^e~^ and intercliange tlie sum- 
mation over m with the integration over z. This yields finally 



9{x) = ,1 i^e-"^'' + ^(1 + erf vW2) (3.40) 

with the error function erf(z) = {2/^) dt e~*^. For x the correlation 
function is dominated by Pi{x) and correspondingly diverges as Xj^Jx^ reflecting 
the tendency of particles to bunch together in the continuous time case. For 
X ^ oo g[x) decays somewhat faster than exponentially, as 

g{x)-p^-J-—e-^-l\ (3.41) 
\J l-Ki^pxy 

Alternatively the correlations between particles can be characterized through 
the variance {ANl)"^ of the number of particles in an interval of size L. When 
L is small compared to the mean interparticle spacing A^^^: is either or 1, and 
(ANlY = pL. For L ^ 1/p a central limit argument shows that 

(ANl^ ^ xL (3.42) 



where the "compressibility" x (defined in analogy with equilibrium systems |]32| ) 
is given by 

xip)=pH{u')~{ur) = p/u, (3.43) 

with the parameter u of the headway distribution ( |1.7[ ). Thus the slope of (ANl)'^ 
versus L changes from unity for L^l/ptol/z/ for L ^ 1/p, reflecting the increase 
(decrease) of particle number fluctuations for continuous time (parallel) dynamics, 
respectively. The compressibility is related to the pair correlation function ( |3.34D 
through 

/oo 
dxigix)-p)). (3.44) 
-oo 
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4 Large scale dynamics of the ARAP 



4.1 Hydrodynamic equation 

The average particle speed v in the ARAP is inversely proportional to the density, 
hence the current j = pv is independent of p. The dynamics on the Euler scale 
X ~ if: is therefore trivial, and one expects a hydrodynamic equation of diffusion 



type |32]. A simple derivation will be given below. Throughout this section we 



consider a general scaled jump length distribution 0(r). 
4.1.1 Continuous time dynamics 

In the continuous time case the ensemble averaged particle positions Xj = (xj) 
evolve according to the linear equations 



dt 



/xi(X,+i-X,). (4.1) 



This problem has been studied previously in the context of crystal growth ||19|| , 
and the procedure can be directly applied to the present context. 

To extract the long wavelength behavior, we introduce a scaling parameter 



3^ , |T5[ e and a smooth function ^(y, r) such that 

X,(t)=e(e2,et). (4.2) 



Inserting this into ( |4.1| ) and expanding to second order in e we obtain 

dr dy 2 dy 



/^i — = — + oTT^- (4-3) 



In the scaling limit e — > this becomes a first order equation which describes 
simple translation to the left [|r^|. 



Here we will however postpone to take the limit, and first carry out a Lagrange 
transformation |2^, |l^, which relates the Lagrangian description in terms of the 
particle positions Xi{t) to the Eulerian evolution of the density field. The local 
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density p near the position of particle i is estimated as (Xj+i — Xj) ^, so using 
we have the relation 



p{ay.r),T) = e-\dildy)-\ (4.4) 

Differentiating this equation with respect to r and using the evolution equation 
]3D for {y, t) one obtains, after some algebra, 

dp dp d f pi\ dp 



e- 



(4.5) 



dt dr dx \2p'^ J dx 
The scaling factor e cancels, and the collective diffusion coefficient is identified to 
be 

Dcip) = 0. (4.6) 

The p~^-dependence is dictated by scale invariance: The typical jump length in a 
region of density p is 6 = pi/ p, and ~ 7(5 ~ p^^. 

4.1.2 Discrete time dynamics 



For discrete parallel update eq.([4.1|) is replaced by 



X,(t+ 1) -X,(t) = /ii[X,+i(t) -X,(t)]. (4.7) 

In the scaling limit e — > this results in the same coarse grained evolution equation 
( [4. 3D , and thus also the nonlinear diffusion equation ([4.5|) is the same as in the 
continuous time case. 

In the case of ordered sequential update one has to take into account that the 
new position of particle z is a random average of its old position and the new 
position of particle i + 1, hence 

X,(t + 1) - X,(t) = /ii[X,+i(t + 1) - X,(t)]. (4.8) 

Making the ansatz Xj(t) = i/p + vt, we see that the average particle speed is 

>^. (4.9) 



p(l - p 
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The speedup compared to continuous time and parallel dynamics is due to the 
decrease of the local density near the update site, see for a discussion of 
similar effects in the asymmetric exclusion process. For the derivation of the 
hydrodynamic equation it is useful to incorporate the expected diffusive scaling 
from the outset and replace ( [4.2| ) by 

X,{t)=aez,eH). (4.10) 



The expansion of ( |4.8| ) to second order in e then yields 



+ (4-11) 



/ii J dr dy 2 dy 

As before, the drift term disappears under the Lagrange transformation based on 
the relation ( [4.4[ ), and one obtains 

^ = = — ( —i^^——] ^ (A 12) 

dt dr dx V2(l -/ii)p2 ) dx' ^ ' ' 

As far as the hydrodynamics is concerned, the different types of dynamics are seen 
to be equivalent up to a rescaling of time. 

4.2 Tracer diffusion 

Hydrodynamic equations of diffusion type are usually associated with symmetric 



(unbiased) particle systems ||3^. In one dimension the tracer diffusion coefficient 
in such systems typically vanishes, and the mean square displacement of a tagged 
particle grows sub diffusively as t^/^ ^. By contrast, the biased random average 
process shows normal tracer diffusion when started from a random initial condition 



and subdiffusive behavior when the initial configuration is ordered [0. Here we 
provide a compact derivation of the two cases and compute the coefficient of the 
asymptotic law for different types of dynamics. 
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4.2.1 Langevin approach for continuous time dynamics 

We start the system in an initial condition without long wavelength fluctuations, 
such as Xi{0) = i/p, i E Z, and denote the positional fluctuation of particle i by 

Ci{t) = Xi{t) - {xi) = Xi{t) - Xi{0) - vt. (4.13) 

For the purpose of extracting the long time behavior of fluctuations, a Langevin 
approximation |jl3[ to the dynamics of Q is suflicient. Thus we add a phenomeno- 
logical noise term rii{t) to the linear equation (|4.1|), 

^ = f^iiQ^i - Q) + m- (4.14) 

The noise is taken Gaussian with zero mean and covariance 

{r],{t)r],{t')) = a6,,6{t-t'). (4.15) 

The noise strength a will eventually be matched to the variance of particle head- 
ways. 

Equation ( [4.14|) is solved by introducing the Fourier transformed fluctuations 



aq,t) = J2e''\n{t) (4.16) 

with wave numbers q in the flrst Brillouin zone [— 7r,7r], and the corresponding 
Fourier transformed noise 

r)(g,t) = ^e^«XW (4-17) 

with covariance 

iviq, mq', t')) = 2naSiq + q')S{t - t'). (4.18) 

The most general quantity of interest is the variance of the displacement between 
particle i at time t and particle j at time t'. By translational invariance this 
depends only on n = i — j and is given by the correlation function 

G„(t,t') = ((Co(t)-Cn(0)')- (4.19) 
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Inserting ( |4.16| ) into ( f4.14| ), solving the equation for C{q,t) and averaging over the 



noise according to ( [4. IS] ) one arrives at the expression 

Gnit, t') = 

1. r ^(2 - - - 2 cos[gn - /x(g)T](e--l^l - e"'^^')) (4.20) 

2vr io ^KQ) 



with uj{q) = /ii(l — cos(g)), /i(g) = /ii sin(g), T = t' — t and T' = t' + t. 

The evaluation is straightforward in the relevant limiting cases. Consider first 
the variance of the headways at time t = t' . For large t ( [4.20|) yields 



Gi(t,t)^-fl-— L=). (4.21) 



III \ 2^nfiit 
This allows us to identify the noise strength a as 



a = fii{{u')-{uf), (4.22) 



and explicitly demonstrates the 1/ v t-approach to the stationary headway distri- 
bution alluded to in ( p.l3| ). 

Next we focus on the dynamics of a single particle and set n = in ( 4.201) . If 



we fix the time increment T = t' — t and let both t and t' oo, Go represents 
the mean square displacement of a particle in the stationary regime. Evaluation 
of ( [4.20| ) gives Go{t,t') ^ cr\T\, which shows that a is precisely the tracer diffusion 



coefficient Dtr- Combining this with ( |4.22| ) and ( |1.5| ) we obtain 

D,, = ^^,{{u') - iu)') = (4.23) 



In fact the first relation in ([4.23|) is easy to understand. The linear equation ([4.3|) 



shows that fiuctuations in the particle positions drift backwards in "label space" 
y = ei. This translates the stationary distance fiuctuations into temporal fiuctu- 
ations, with a conversion factor given by the drift speed fii. As was mentioned 
already, the existence of a nonvanishing tracer diffusion coefficient for models with 
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a hydrodynamic equation of diffusion type is unusual in one dimension, since gener- 
ically such an equation implies symmetric particle jumps, in which case the tracer 
particle displacement grows only sub diffusively due to the single file constraint 
P, 1^. Here Z^tr is nonzero because the particles move, at speed v, relative to the 
(stationary) density fluctuations. A rigorous derivation of ( [4.23| ) has recently been 
presented by Schiitz p9| . 

Since the hydrodynamic equations in the two cases are identical, the argument 
leading to first relation in ([4. 231) carries over directly to discrete parallel update. 



and using (|3.22|) we conclude that the tracer diffusion coefficient in this case is 
given by 

or = ^^ii^^^. (4.24) 



Similarly the expression 



rjseq _ /il(/i2 -/ij) . . 



is obtained for the backward sequential case by combining eqs.( p.29"D and ( [4.9| ). 
Both (WM) and ( ^^25]) have been verified numerically for the case of uniform (f){r). 



Subdiffusive behavior is found in the mean square displacement of a particle 
starting from an initial configuration without long wavelength disorder []T^. This 
is given by (^4.20|) with n = t' = 0. For large t one obtains 



i,m^rM,0)..,l±-^-,^^^'-^. (4.26) 



Using ( |3.43| ) and ( [4.6| ) this is seen to agree with the expression 

{cm = vv^ix/p')VD:t (4.27) 

derived from hydrodynamic arguments 

4.2.2 The independent jump approximation 

For the totally asymmetric simple exclusion process it is known [^, ^ |ll[ that the 
motion of a tagged particle in the stationary state follows a Poisson process, and 
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therefore the tracer diffusion coefficient is simply equal to the mean speed 1 — p. 
Here we show that the expressions ( [4.23| - [4. 251 ) for the ARAP are consistent with 



a similar independent jump picture. 

Consider first the case of discrete time dynamics, where the random choice 
of the jump length 5j is the only source of disorder, and therefore the tracer 
diffusion coefficient for independent jumps is equal to the variance of 5j. For 
parallel update 5i is a uniform random fraction of the particle headway Wj, hence 
(5^) — (5)^ = /i2(^i^) — Pi/p^j which is easily checked to coincide with ( 4.24| ). For 



the backward sequential case 5i is a random fraction of the intermediate state 
headway u[. Therefore, using eqs. ( |3.28| ), ( |3.27| ) and ( p.29| ). 



which is also found to agree with ( [4.25|) . 



f). (4-28) 



In the continuous time case the random timing of jumps introduces an addi- 
tional source of disorder. It is natural to assume, in analogy with the asymmetric 
exclusion process, that the jumps occur according to a Poisson process. In the 
independent jump approximation the particle displacement Ax in time t is then 
given by 

n(t) 

Aa;(t) = ^ <5(') (4.29) 
1=1 

where n{t) is a Poisson random variable with mean t and the jump lengths 6^^^ are 
independent random fractions of the (independent, random) particle headways. It 
is straightforward to show that the variance of Ax is 

((Ax)2) - {Axf = {6^)t, (4.30) 

thus in this case the independent jump approximation to D^r is (5^) = ^2{u^) in 
agreement with (|4.23| ). 
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5 Summary and outlook 



We have presented results for two classes of particle systems on M. The models 
considered in Section 2. have Poisson invariant measures and nonlinear current- 
density relations (see eqs.( p.8| , p.l4| )). Time-dependent fluctuations in these models 



are therefore expected to be governed by the noisy Burgers (or Kardar- 

Parisi-Zhang |]14|) equation, which is not amenable to simple analysis. By contrast, 



the asymmetric random average processes introduced in Section 3. have nontrivial 
invariant measures, but the linearity of the jump rules allows for a detailed study 
of dynamic properties (Section 4.). 

A central result for the ARAP is the dependence of the headway distribution 
( |L7|) on the type of dynamics. The idea that parallel update reduces density 
fluctuations is familiar from earlier work on the asymmetric exclusion process and 
related models for traffic flow, however in that case the ordered sequential update 
produces the same (Bernoulli) invariant measure as the continuous time process 



Our study suggests that the invariant measure of the continuous time ARAP 
displays an unusual combination of features: The two-point headway correlations 
factorize, the single particle headway distribution appears to be exactly given by 
the expression ( |1.7| ) derived under the assumption of pairwise independence, but 
nevertheless the product measure (|1.3| ) is not invariant. Rajesh and Majumdar 
have found the same features in a larger class of models which interpolate between 
continuous time and parallel update It would be most interesting to find a 



simple "deformation" of the product measure which explains this behavior. The 
status of the product measure assumption for the ordered sequential update also 
remains to be clarified. The considerations of Section 3.2.2 indicate that it might 
be possible to exactly reduce this case to that of parallel update, for which the 
product measure is known to be invariant 0. 
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Another interesting direction for future work is the introduction of quenched 
random inhomogeneities. In asymmetric exclusion models it is possible to find 
invariant product measures also in the presence of random jump rates associated 
with particles p, ||, For the continuous time ARAP with jump rates 7^ 
depending on the particle label i (the position i in the stick representation) pre- 
liminary numerical simulations indicate that the product measures discussed above 
do not persist. It is possible to write down a closed set of linear equations for the 
two-point function (uiUj) which depends on the disorder configuration {74} and 
which should yield insight into the emergence and nature of correlations. Here 
we merely remark that, since the mean speed of particle i is 7i(Mj), stationarity 
implies (ui) = C/'-fi where the constant C is fixed by the average headway. If 
the distribution of jump rates is chosen such that (l/7i) exists, C (l/'-fi)^'^ in 
the limit of infinite system size, and all headways have a finite mean. Otherwise 
(e.g. for a uniform distribution of jump rates) arbitrarily large headways will open 
in front of the slowest particles, similar to the low density phase of asymmetric 
exclusion models with particlewise disorder |21|, ||, |^ . 
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Dynamics: 






\ 1 




Continuous 


2.998 ± 1 
(3) 


15.02 ± 1 

(15) 


105.3 ±2 
(105) 


947 ±3 
(945) 


Sequential 


1.9997 ± 2 

(2) 


5.995 ±2 
(6) 


23.93 ±2 
(24) 


119.2 ±2 
(120) 


Parallel 


1.4998 ± 1 

(3/2) 


2.996 ± 1 

(3) 


7.466 ± 3 
(15/2) 


22.26 ±2 

(45/2) 



Table I. The Table contains numerical estimates of the first few moments of the 
stationary headway distribution for the ARAP with uniform 0(r) and different 
kinds of update. The data were obtained from simulations of systems of 2 x 10^ 
particles which were started from an ordered initial condition, = 1 for all i, 
and allowed to evolve for lO'' time steps. To extrapolate to t — > oo, each run was 
fitted to eq. (|3.13 ), and the errors were estimated by taking an average over 10 
runs (errors refer to the last digit shown). The numbers in parentheses are the 
conjectured values of the moments; for the case of parallel update these are known 
to be exact 0. The remaining discrepancies are in fact largest for parallel update, 
and can probably be attributed to residual finite time effects. 
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